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The confirmatory analysis of pre-specified multiple hypotheses has become common in pivotal 
clinical trials. In the recent past multiple test procedures have been developed that reflect the relative 
importance of different study objectives, such as fixed sequence, fallback, and gatekeeping proce- 
dures. In addition, graphical approaches have been proposed that facilitate the visualization and 
communication of Bonferroni-based closed test procedures for common multiple test problems, such 
as comparing several treatments with a control, assessing the benefit of a new drug for more than one 
endpoint, combined non-inferiority and superiority testing, or testing a treatment at different dose 
levels in an overall and a subpopulation. In this paper, we focus on extended graphical approaches by 
dissociating the underlying weighting strategy from the employed test procedure. This allows one to 
first derive suitable weighting strategies that reflect the given study objectives and subsequently apply 
appropriate test procedures, such as weighted Bonferroni tests, weighted parametric tests accounting 
for the correlation between the test statistics, or weighted Simes tests. We illustrate the extended 
graphical approaches with several examples. In addition, we describe briefly the gMCP package in R, 
which implements some of the methods described in this paper. 

Keywords: Dunnett test; Gatekeeping procedure; Min-p test; Non-inferiority; Truncated 
Holm. 



1 Introduction 

Multiple test procedures are often used in the analysis of clinical trials addressing multiple objec- 
tives, such as comparing several treatments with a control and assessing the benefit of a new drug 
for more than one endpoint. Several multiple test procedures have been developed in the recent past 
that allow one to map the relative importance of the dilferent study objectives as well as their 
relation onto an appropriately tailored multiple test procedure. 

A common strategy to reduce the degree of multiplicity is to group the hypotheses into primary 
and secondary objectives (O'Neill, 1997). Test procedures accounting for the inherent logical 
relationships include fixed sequence tests (Maurer et al., 1995; Westfall and Krishen, 2001), gate- 
keeping procedures (Bauer et al., 1998; Westfall and Krishen, 2001; Dmitrienko et al., 2003) and 
fallback procedures (Wiens, 2003; Huque and Alosh, 2008). Li and Mehrotra (2008) introduced a 
more general approach for adapting the significance level to test secondary hypotheses based on the 
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finding for tiie primary iiypotlieses. Alosli and Huque (2009) introduced tlie notion of consistency 
when testing for an effect in the overall population and in a specific subgroup. The authors extended 
this consistency concept to other situations (Alosh and Huque, 2010), including how to address 
multiplicity issues of a composite endpoint and its components in clinical trials (Huque et al., 2011). 
Hung and Wang (2009, 2010) considered some controversial multiple test problems, with emphasis 
on regulatory applications, and pointed out illogical problems that may arise with recently devel- 
oped multiple test procedures. 

In this paper, we focus on graphical approaches which have been introduced independently by 
Bretz et al. (2009) and Burman et al. (2009). The key idea is to express the resulting multiple test 
procedures by directed, weighted graphs, where each node corresponds to an elementary hypothesis, 
together with a simple algorithm to generate such graphs while sequentially testing the individual 
hypotheses. Using graphical approaches, one can explore different test strategies together with the 
clinical team and thus tailor the multiple test procedure to the given study objectives. So far, the 
description of these graphical approaches has focused on Bonferroni-based test procedures. In this 
paper, we investigate extensions of the original ideas. In particular, we discuss in Section 2 how a 
separation between the weighting strategy and the test procedure facihtates the application of a 
graphical approach beyond Bonferroni-based test procedures. In Section 3, we illustrate these ideas 
with different test procedures. We start with a brief review of Bonferroni-based test procedures and 
subsequently describe parametric graphical approaches that account for the correlation between the 
test statistics as well as graphical approaches using the Simes test. In Section 4, we describe the 
gMCP package in R which implements some of the methods discussed in this paper and illustrate it 
with a clinical trial example using a truncated Holm procedure. Concluding remarks are given in 
Section 5. 



2 Graphical weighting strategies 

Consider the problem of testing m elementary hypotheses i/i,...,//^, some of which could be more 
important than others, e.g. primary and secondary objectives. Let /={!,.. .,m} denote the asso- 
ciated index set. The closure principle introduced by Marcus et al. (1976) is commonly 
used to construct powerful multiple test procedures. Accordingly, we consider all non-empty 
intersection hypotheses Hj — D/^jIIj, J c. I. We further pre-specify an a-level test for each 
Hj. The resulting closed test procedure rejects Hi, i e /, if all intersection hypotheses Hj with 
i & J c. I are rejected by their corresponding a-level tests. By construction, closed test procedures 
control the familywise error rate (FWER) in the strong sense at level ae(0,l). That is, the prob- 
abihty to reject at least one true null hypothesis is bounded by a under any configuration 
of true and false null hypotheses (Hochberg and Tamhane, 1987). In fact, closed test procedures 
have certain optimality properties whenever the FWER has to be controlled (Bauer, 1991). In what 
follows, we assume that the hypotheses H\,...,Hm satisfy the free combination condition (Holm, 
1979). If this condition is not satisfied, the methods in this paper still control the FWER at level a, 
although they can possibly be improved because of the reduced closure tree (Brannath and Bretz, 



One important class of closed test procedures is obtained by applying weighted Bonferroni tests 
to each intersection hypothesis Hj. For each J^I assume a collection of weights Wj{J) such that 
0 < Wj{J) < 1 and ^ 1- With the weighted Bonferroni test we reject Hj if pj < olj(J) = 

Wj(J)a for at least one /e/, where pj denotes the unadjusted p-value for Hj. Hommel et al. (2007) 
introduced a useful subclass of sequentially rejective Bonferroni-based closed test procedures. They 
showed that the monotonicity condition 



2010). 



Wj(J) < Wj(J') for all / c / c / and j e /' 



(1) 
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ensures consonance, i.e. if an intersection hypothesis Hj is rejected, there is an index jeJ, such that 
the elementary hypothesis Hj can be rejected as well. This substantially simplifies the 
implementation and interpretation of related closed test procedures, as the closure tree of 2™—! 
intersection hypotheses is tested in only m steps. Many common multiple test procedures satisfy (1), 
see Hommel et al. (2007) for examples. 

Bretz et al. (2009) and Burman et al. (2009) independently derived graphical representations and 
associated rejection algorithms for important subclasses of the Hommel et al. (2007) procedures. The 
graphical representations and rejection algorithms in these two articles are different, though under- 
lying ideas are closely related; see Guilbaud and Karlsson (2011) for some comparative examples. 
Using the graphical approach of Bretz et al. (2009), the hypotheses H H„ are represented by 
vertices with associated weights denoting the local significance levels ai,...,a„,. In addition, any two 
vertices if, and Hj are connected through directed edges, where the associated weight gy indicates the 
fraction of the (local) significance level a,- that is propagated to Hj once Hi (the hypothesis at the tail of 
the edge) has been rejected. A weight gij = 0 indicates that no propagation of the significance level is 
foreseen and the edge is dropped for convenience. Figure 1 shows an example. 

While the original graphical approaches were introduced based on weighted Bonferroni tests, we 
propose here to dissociate the underlying weighting strategy from the employed test procedure. The 
benefit of such an approach is the enhanced transparency by (i) first deriving suitable weighting 
strategies that reflect the given study objectives (and which can be communicated to the chnical 
team) and (ii) subsequently applying appropriate test procedures that do not necessarily have to be 
based on Bonferroni's inequality. 

Graphical weighting strategies are conceptually similar to the graphs proposed by Bretz et al. 
(2009). They essentially summarize the complete set of Ym=\ ' (D = ^2™"' weights determining the 
full closure tree. A weighted multiple test can then be apphed to each intersection hypothesis Hj, 
such as a weighted Bonferroni test, a weighted min-p test accounting for the correlation between the 
test statistics, or a weighted Simes test; see Section 3 for details. Weighting strategies are formally 
defined through the weights iv,(/), iel, for the global null hypothesis Hj and the transition matrix 
G = {gij}, where 0<^,y< 1, gu = 0, and Yl"j'=\ Sn < 1 for all We additionally need to determine 

how the graph is updated once a vertex is removed. This can be achieved by tailoring Algorithm 1 in 
Bretz et al. (2009) to the graphical weighting strategies as foUows. For a given index set /^/, let 
J" = I\J denote the set of indices that are not contained in /. Then the foUowing algorithm de- 
termines the weights \Vj{J), jeJ. This algorithm has to be repeated for each J^I to generate the 
m2'"~^ weights for the full closure. 

Algorithm 1 (Weighting Strategy) 

(i) Select jeJ^ and remove Hj 

(ii) Update the graph: 



(iii) If |/'| > 1, go to step (i); otherwise W((,J) — wt{I),i e /, and stop. 

As shown by Bretz et al. (2009), the weights Wj{J), jeJ are unique. In particular, they do not 
depend on the sequence in which hypotheses Hj, j e are removed in step (i) of Algorithm 1. Note 
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Figure 1 Weighting strategy for two hierarchically ordered endpoints and two dose 
levels. 

that Algorithm 1 requires specifying the weights 117(7) for the global intersection hypothesis Hj and 
the elements of the transition matrix G. This leads to the specification of m+m(m — \) — tr? 
parameters if X^ye/ ^'A^) — 1 ^^'^ Y^jLi Sa < 1 or — l+m{m — 2) — m^ — m—\ parameters if 
H'//) = 1 and = 1, for all Ujel. 

Example 1 

As an example, assume a primary family of two hypotheses — {7/1,7/2} and a secondary 
family of two hypotheses J^i — {7/3,7/4}. The hypotheses 7/i and H2 could denote, for example, the 
comparison of low and high dose with a control, for either a primary endpoint, a non-inferiority 
claim, or an overall population. Accordingly, the hypotheses 7/3 and 7/4 would then denote the 
comparison of the same two doses with a control, for either a secondary endpoint, a superiority 
claim, or a pre-specified subgroup. Figure 1 visualizes one possible weighting strategy. It is moti- 
vated by a strict hierarchy within dose: the secondary endpoint will only be assessed if efficacy was 
shown previously for the primary endpoint (so-called successiveness property; see Maurer et al., 
2011). If for one of the doses efficacy can be shown for both the primary and the secondary 
endpoint, the associated weight is passed on to the other dose. Therefore we have I — {1,2,3,4}, 
vi'i(/) = W2{I) — 0.5 for the primary hypotheses and 1^3(7) — 11*4 (/) — 0 for the secondary hy- 
potheses, which imphes that no secondary hypothesis can be rejected until a primary hypothesis is 
rejected and propagates its weight. The associated transition matrix is 

/O 0 1 0\ 

0 0 0 1 

0 10 0" 
yi 0 0 0/ 

The graph in Figure 1 together with Algorithm 1 from above fully specify the 32 weights of the 
full closure tree, as summarized in Table 1. This table parallels the weight tables introduced by 
Dmitrienko et al. (2003). Note that the weights Hy(/), j e are formally not defined and expressed 
by "-" in Table 1. Figure 2 displays the updated graphs resulting from Figure 1 after removing Hi, 
7/2, /73, or 7/4. The four updated graphs in Figure 2 correspond to the four rows in Table 1 
containing the weights for the three-way intersection hypotheses. Removing any two hypotheses 
results in six possible two-way intersection hypotheses and the two vertexes are connected by two 
directed edges, each with weight 1 (graphical display omitted here). Note that Figure 2 displays the 
principle of recalculating the weights by updating the graphs. It is possible and also necessary to 
remove hypotheses with weight 0 (in this example 7/3 and H4 with h'3(7) — m'4(/) — 0) in order to 
compute the respective weights for the larger intersection hypotheses. 

Note that Figure 1 displays only one possible weighting strategy. Many other weighting strategies 
are possible and perhaps more reasonable, depending on the given context. We refer to Bretz et al. 
(2011) for a generic discussion about testing two families J^i and ^2 with two hypotheses each. 
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Table 1 Weights for the intersection hypotheses derived from Figure 1. 



Intersection hypothesis Weights 





Hi 


H2 


H3 


H4 


HinH2nH2nH4 


0.5 


0.5 


0 


0 




0.5 


0.5 


0 






U.J 


U.J 




A 
U 




0.5 


0.5 






HinHj,nH4 


0.5 




0 


0.5 




1 




0 




//,n//4 


0.5 






0.5 


Hi 


1 








H2nHjr\H4 




0.5 


0.5 


0 






0.5 


0.5 




H2nH4 




1 




0 


H2 




1 






//,n//4 






0.5 


0.5 








1 




H4 








1 




Figure 2 Updated graphs resulting from Figure 1 after removing (A) Hi, (B) H2, 
(C) H3, and (D) H4. 

3 Test procedures 

In Section 2, we proposed to dissociate the underlying weighting strategy from the employed test 
procedure and gave a generic description of the former, illustrated with an example. In this section 
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we give details on dilferent test procedures that could be employed to test the intersection 
hypotheses, including weighted Bonferroni tests, weighted min-;? tests accounting for the correlation 
between the test statistics, and weighted Simes' tests. 



3.1 Weighted Bonferroni tests 

The weighted Bonferroni test introduced in Section 2 is the simplest apphcable test procedure, 
leading to the original graphical approaches by Bretz et al. (2009). Applying the Bonferroni test 
leads to simple and transparent test procedures that are often easier to communicate than alter- 
native, potentially more powerful approaches. As a matter of fact, the Bonferroni test is often 
perceived to provide credible trial outcomes in chnical practice. Most importantly in the context of 
the graphical weighting strategies considered here, applying the Bonferroni test leads to shortcut 
procedures as long as the monotonicity condition (1) is satisfied. That is, one can start with a graph 
as shown in Figure 1 and sequentially test the m hypotheses as long as individual null hypotheses //,-, 
/ el, are rejected. Based on Algorithm 1 from Section 2, we give in the following a similar algorithm 
that accounts for the weighted Bonferroni tests, thus leading to the sequentiaUy rejective multiple 
test procedures described in Bretz et al. (2009): 



Algorithm 2 (Weighted Bonferroni Test) 

(i) Select a jel such that Pj<Wj{I)a. and reject H/, otherwise stop. 

(ii) Update the graph: 



/ ^ AW 

r weil) + wjil)gje, {€l 
\ 0, otherwise 



W(il) 

ga ^ 



0, otherwise 



(iii) If |/| > 1, go to step (i); otherwise stop. 

Similar to Algorithm 1, the results in Bretz et al. (2009) ensure that the decisions of the resulting 

sequentially rejective multiple test procedures remain unchanged regardless of the actual rejection 
sequence. That is, if in step (i) of Algorithm 2 more than one hypothesis could be rejected, it does 
not matter with which to proceed. Although Algorithms 1 and 2 have a similar update rule in 
step (ii), they differ in the way that the index sets are updated. While Algorithm 2 starts with the 
global index set / and reduces it sequentially as long as hypotheses are rejected, Algorithm 1 
removes, for each ./^/, consecutively all indices from / that are not contained in / until the set J is 
obtained. Note that performing a closed weighted Bonferroni test procedure using the weights from 
Algorithm 1 leads to exactly the same test decisions as performing a sequentially rejective multiple 
test procedure with Algorithm 2 based on the same starting weights. 

Figure 3 gives an example of a Bonferroni-based sequentially rejective multiple test procedures 
for the weighting strategy proposed in Example 1. Assume, for example, the unadjusted /7-values 
Pi = 0.01, p2 = 0.005, pj, = 0.1, and p4 = 0.5. Then we can reject both Hi and Hj, but none of the 
other hypotheses. Figure 3 displays the initial graph together with a possible rejection sequence. As 
mentioned above, the final decisions on which hypotheses to reject do not depend on the particular 
rejection sequence. That is, with the initial graph from Figure 3 we would obtain the same decisions, 
regardless of whether we first reject H2 and then Hi, or vice versa. 
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initial graph 



H2 rejected 



H2, Hi rejected 



OL 



a 



2 



2 



2 




Figure 3 Graph for sequentially rejective procedure with example rejection sequence. 

Many standard approaches from the hterature can be visualized using Bonferroni-based 
graphical test procedures, including the weighted or unweighted Bonferroni-Holm procedure 
(Holm, 1979), fixed sequence tests (Maurer et. al, 1995; Westfall and Krishen, 2001), fallback 
procedures (Wiens, 2003), and gatekeeping procedures (Bauer et al., 1998; Westfall and Krishen, 
2001; Dmitrienko et al., 2003). Adjusted ;7-values and simultaneous confidence intervals can be 
calculated as well, although the resulting simultaneous confidence intervals are known to be of 
limited practical use, as they are often non-informative; see Strassburger and Bretz (2008), Guilbaud 
(2008, 2009) and Bretz et al. (2009) for details. Bretz et al. (2011) provided SAS/IML code to 
perform the resulting Bonferroni-based sequentially rejective multiple test procedures. In Section 4, 
we describe the gMCP package in R, which offers a convenient graphical user interface (GUI) for 
these approaches. 

One general disadvantage of Bonferroni-based approaches is a perceived power loss, motivating 
the use of weighted parametric tests that account for the correlation between the test statistics or the 
use of weighted Simes tests. We discuss these alternative test procedures in Sections 3.2 and 3.3, 
respectively. 

3.2 Weighted parametric tests 

If for the intersection hypotheses Hj,J c /, the joint distribution of the p-values pj, jeJ, are known, 
a weighted min-p test can be defined (Westfall and Young, 1993; WestfaU et al., 1998). This test 
rejects Hj if there exists a jeJ such that pj < cjWj{J)a, where Cj is the largest constant satisfying 



If the ;7-values are continuously distributed, there is a cy such that the rejection probability is 
exactly a. Determination of cy requires knowledge of the joint null distribution of the ;;-values and 
computation of the corresponding multivariate cumulative distribution functions. If the test sta- 
tistics are multivariate normal or t distributed under the nuU hypotheses, these probabihties can be 
calculated using, for example, the mvtnorm package in R (Genz and Bretz, 2009). Alternatively, 
resampling-based methods may be used to approximate the joint null distribution; see Westfall and 
Young (1993). 

If Cy= 1 in (2), the weighted parametric test reduces to the weighted Bonferroni test. This fully 
exhausts the level if and only if the joint distribution of continuously distributed ;7-values with 
strictly positive density function over (0,1)"' satisfies 




(2) 



PHj(.{Pj < cjWjiJ)a} n {pt < cjWiiJ)a}} = 0 
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for all i^jeJ, because then all events are pairwise disjoint and Puji^jeJ [pj < cjWj{J)'x}) 
= Phj(Pj < cjWj(J)ol). Otherwise, Cj>l and the weighted parametric test gives a uniform im- 
provement over the weighted Bonferroni test from Section 3.1. 

If not all, but some of the multivariate distributions of the ;7-values are known, it is possible to 
derive conservative upper bounds of the rejection probability that still give an improvement over the 
Bonferroni test. Assume that / can be partitioned into / sets //, such that / — ^'/j^ilh and /,■ n //, = 0 
for i ^ h— 1 , . . ., /. We assume that for each h= 1 , . . . ,/ the joint distribution of the />-values pi, i e Ih, 
is known, but the joint distribution of /^-values belonging to different 1^ is not necessarily known. 
Now, let J^I and choose the maximal critical value cj such that 

Y.PMV\^k< c/Wfc(y)a} ) < a. (3) 

h=\ \fre/,,n./ / 

By the Bonferroni inequahty, the left-hand side in (2), which cannot be computed if the full joint 
distribution is unknown, is bounded from above by the left-hand side in (3), whose computation 
requires only the knowledge of the joint distribution of the /7-values in //,n/, separately for each 
h= 1,...,/. Thus, any cj satisfying (3) will also satisfy (2), leading to a conservative test for the 
intersection hypothesis Hj. 

It follows immediately from Eq. (1) that these parametric approaches are consonant if 

cjWj(J) < cj'Wj(J') for all / c / c / and ./ e f. (4) 

For values following a joint continuous distribution with strictly positive density function 
over (0,1)'" this is also a necessary consonance condition. This condition is often violated by the 

weighted parametric tests above. Consider, for example, the Sidak (1967) test for three hypotheses 
with initial weights 1/3. Assume that for the test of the intersection of any two hypotheses the 
weights are 1/3 and 2/3. For a = 0.05, the critical value CjWj (/)oc = 0.01695 for all three hypotheses 
in the first step. For all / with |/| = 2, we have cy'H',(/')a = 0.01686 for the hypothesis with the 
weight 1/3 in the second step, violating (4). This phenomenon is even more pronounced for positive 
correlations. If in the previous example the correlations are all 0.5 (corresponding to a Dunnett test 
in a balanced one-way layout with known variance), we have cjWj(J)a. = 0.0196 and 
crWjiJ')a = 0.0182. 

If the consonance condition (4) is met, a sequentially rejective test procedure similar to the 
Bonferroni-based graphical tests from Section 3.1 can be defined. 

Algorithm 3 (Weighted Parametric Test) 

(i) Choose the maximal constant cj that satisfies either (2) or (3) for / = /. 

(ii) Select a jel such that pj < ciWj{I)a. and reject Hf, otherwise stop. 

(iii) Update the graph: 



gek 



A{/} 

^ r we{I)+Wj{I)gje, (€l 
\ 0, otherwise 

\ 0, otherwise 



(iv) If |/| > 1, go to step (i); otherwise stop. 

For any specific multiple test procedure defined by a given graph, the consonance condition can 
be checked. If the consonance condition is not met, the weighting strategies introduced in Section 2 
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remain applicable, although the connection to a corresponding sequentially rejective test procedure 
is lost. In this case, Algorithm 3 no longer appUes and one has to go through the entire closed test 
procedure. For a given weighting strategy, this procedure is uniformly more powerful than the 
associated Bonferroni-based procedure from Section 3.1. Note that adjusted /^-values for each 
hypothesis //, can be obtained by computing ;7-values for each intersection hypothesis Hj with ieJ 
(given by the lowest local level for which the respective intersection hypothesis can be rejected) and 
then taking the maximum over them. 

Before illustrating Algorithm 3 with two examples, we notice that Eq. (2) does not provide the 
only possible definition of a weighted parametric test. Instead of using cjWj{J)a as the critical values 
for pj, je J, we could also use some other function //(vt',(/),a) fulfilling /y(i47(/), a) > Wj{J)a for all 
jeJ and all dependence structures of the /^-values. For example, if Tj — $"'(1 — pj) is the test 
statistic corresponding to the /?- value of a z-test for Hj, then finding an 8/ such that 



would also define a test which is uniformly more powerful than the corresponding weighted Bon- 
ferroni test. A related approach to account for correlations in weighted multiple testing procedures 
defined by the graphical approach was considered in Millen and Dmitrienko (2011). 



We revisit the weighting strategy from Example 1. Assume that the joint null distribution 
of the /^-values pi, p2 for the two primary dose-control comparisons as well as the joint null 
distribution of the /(-values pj„ p^ for the two secondary comparisons are known. Applying the 
standard analysis-of-variance assumptions with a known common variance, we have a bivariate 
normal distribution, where the correlation is determined only by the relative group sample sizes. In 
practice, the correlation between primary and secondary endpoints is typically unknown and thus 
the joint distributions of the pairs (pi,Pj),i = 1,2, j = 3,4 are also unknown. Therefore, (2) cannot 
be computed and C/ cannot be determined directly. Setting I\ — {1,2} and I2 — {3,4}, the joint null 
distribution of the test statistics for the hypotheses in Ii and I2 is known and the constants cj can be 
determined by (3). Note that cj depends on a and on the weights. Table 2 shows the local sig- 
nificance levels for both (A) the closed weighted Bonferroni test procedure and (B) the closed 
weighted parametric test procedure, assuming a = 0.025 and equal group sample sizes. 

Using, for example, the mvtnorm package in R, one can call 

> myf ct <- function(x, a, w, sig) { 

-I- 1 - a - pmvnorm ( lower = -Inf , upper = qnorm ( l-x*w*a) , sigma = sig) 



> sig <- diag{2)*0.5 -1-0.5 

> uniroot (myf ct , lower = 1 , upper = 9,a = 0.02 5,w = rep (0.5, 2 ) , 
+ sig = sig) $root 

[1] 1.078306 

to compute C/= 1.0783 for J — {3,4} as well as for all /^{1,2} and Cj=l otherwise. 
In other words, H2,f\H^ and all intersection hypotheses that include Hi and H2 are 
tested with unweighted Dunnett z tests. However, intersection hypotheses containing H\r\Hi or 
H2C^H^ are tested with an unweighted Bonferroni test. As a consequence, the resulting family of 
tests is not consonant. For example, C(|^2,3,4)M'i({1, 2,3, 4})a = 0.0135>0. 0125 = Cji,4)H'i({l,4})a, 
violating condition (4). Nevertheless, for a given weighting strategy, the closed test 
procedure based on parametric weighted tests dominates the associated procedure based on 
weighted Bonferroni tests. For example, if /ii = 0.0131, ;?2 = 0.1, = 0.012, and p4 = Q.0\, the 




Example 2 



+ } 
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Table 2 Local significance levels (in %) of A: weighted Bonferroni (B: parametric, C: consonant 
parametric with 5 = 0.0783) test for the example from Figure 1 and a = 0.025. 



Intersection hypothesis Local significance levels (in %) 





Hi 


H2 


Hs 


H4 


TJ r^TJ r^Tl r\TJ 
rlyl 1/22' \-tl4 


l.Zj (i.3j,i.Jjj 


l.Zj (i.JJ,i.3j) 


A /A A\ 

0 (0,0) 


A /A A\ 

0 (0,0) 


HinH2nH3 


1.25 (1.35,1.35) 


1.25 (1.35,1.35) 


0 (0,0) 




Hir\H2r\H4 


1.25 (1.35,1.35) 


1.25 (1.35,1.35) 




0 (0,0) 


Hir\H2 


1.25 (1.35,1.35) 


1.25 (1.35,1.35) 






Hir\H3nH4 


1.25 (1.25,1.35) 




0 (0,0) 


1.25 (1.25,1.15) 


HiDHs 


2.50 (2.50,2.50) 




0 (0,0) 




HinH4 


1.25 (1.25,1.35) 






1.25 (1.25,1.15) 


Hi 


2.50 (2.50,2.50) 












1.25 (1.25,1.35) 


1.25 (1.25,1.15) 


0 (0,0) 


7/207/3 




1.25 (1.25,1.35) 


1.25 (1.25,1.15) 




7/207/4 




2.50 (2.50,2.50) 




0 (0,0) 


H2 




2.50 (2.50,2.50) 






7/307/4 






1.25 (1.35,1.35) 


1.25 (1.35,1.35) 


7/3 






2.50 (2.50,2.50) 




7/4 








2.50 (2.50,2.50) 



weighted parametric test procedure rejects 7/i and 7/3, whereas the Bonferroni test rejects none. In 

Section 4, we revisit this numerical example and describe the gMCP package in R, which implements 
the closed weighted parametric test procedure (B). Related gatekeeping procedures addressing the 
problem of comparing several doses with a control for multiple hierarchical endpoints were 
described, among others, by Dmitrienko et al. (2006), Liu and Hsu (2009), and Xu et al. (2009). 

Continuing with the example, one can enforce consonance via an appropriate modification 
of the weighting strategy from Figure 1. To achieve consonance, we introduce additional 
edges with weight 5 (see Figure 4) such that the weight for Hi (resp. 7/2) is sufiiciently increased to 
satisfy the mono tonicity condition (4) when testing the intersection hypotheses 7/1O7/4 and 
7/1O7/3O7/4 (resp. 7/207/3 and //2O//3O7/4). If § > 8*: — Cji,2,3,4| - 1 the resulting closed 
test procedure is consonant and Algorithm 3 can be used to perform the test. In the above 
example with a = 0.025, the lower bound is 5* = 0.0783. Setting 5 = 5*, we obtain the local sig- 
nificance levels for procedure (C) in Table 2. Note that because of the special weighting strategy 
employed in this example, these local significance levels are obtained with the regular Dunnett and 
univariate z tests. 

The lower bound S* depends on the correlation between the test statistics for Hi and //2. Because 
C{i,2,3,4) increases with the correlation, this also holds for 8*. In the limiting case that the sample 
size ratios of the dose groups and the control group tend to infinity, the correlation tends 
to 1. Consequently, C(i_2,3,4) = 2, such that 5* = 1 and the graph is degenerated for all a>0. On the 
other hand, if the above sample size ratios tend to 0, the correlation tends to 0 and 5* = 
2(1 -(1 -a)'/2)/a- 1 in limit. 

Note that by enforcing consonance, the resulting multiple test procedure based on weighted parametric 
tests is no longer uniformly better than the associated Bonferroni-based test procedure which does not 
account for the correlations. That is, for a given weighting strategy, the closed test procedure based on 
parametric weighted tests may fail to reject certain hypotheses that otherwise are rejected by the asso- 
ciated procedure based on weighted Bonferroni tests. For example, if pi — 0.01, ;72 = 0.1, 773 = 0.012, 
and /74 = 0.01, the initial graph from Figure 3 rejects Hi and H^, whereas the consonant weighted 
parametric test procedure from Figure 4 with 5 = 0.0783 rejects only Hi. 
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Figure 4 Graphical display of weighting strategy for a consonant weighted parametric 
test procedure. 

Example 3 

Consider again Example 1, but assume that Hx^Hj are two non-inferiority hypotheses (say, for 
low and high dose against control) and H3,H4 are two superiority hypotheses (for the same two 
doses). We again make the standard analysis-of-variance assumptions with a known common 
variance and let a = 0.025. Bonferroni-based graphical approaches for combined non-inferiority 
and superiority testing were described in Hung and Wang (2010) and Lawrence (2011). In the 
following, we exploit the fact that the correlations between the four test statistics are known. 
Therefore, the complete joint distribution is known and we can apply (2). Note that if uy(/) = 0 for 
some /e/, the joint distribution degenerates. In our example it thus suffices to calculate bivariate or 
univariate normal probabihties. 

Assume first that the same population is used for all four tests. For simplicity, assume further that 
the group sample sizes are equal. Then the correlation between the non-inferiority and superiority 
tests within a same dose is 1; all other correlations are 0.5. Therefore, Cj= 1.0783 for 
/ = {1,2}, {1,4}, {2, 3}, and {3,4}. Otherwise, cj= 1 and condition (4) is trivially satisfied. That is, 
consonance is ensured and one can apply Algorithm 3. This leads to a sequentially rejective multiple 
test procedure, where at each step either bivariate Dunnett z tests or individual z tests are used. This 
conclusion remains true if the common variance is unknown and Dunnett t tests or individual t tests 
are used. 

To illustrate the procedure, let a = 0.025 and assume the unadjusted ;;-values pi = 0.01, />2 = 0.02, 
P3 = 0.005, and p4 = 0.5. Following Algorithm 3, we have pi < c/M'i(/)a — 0.0135 and can reject Hi. 
The update step then leads to the weights in Figure 2(A). Next, p^ < 0.0135 and we can reject H^. 
This leaves us with H2, H4 and the weights h'2({2, 4}) = 1, h'4({2, 4}) — 0. Therefore, Hj is now tested 
at full level a. Because P2<^, we reject H2 and the procedure stops. 

We now consider the situation that two different populations are used. Assume that the per- 
protocol population (PP) is used for non-inferiority testing and the intention-to-treat population 
(ITT) for superiority testing, where PP is a subpopulation of ITT. Let «, denote the ITT sample size 
for group where / = 0 (1,2) denotes placebo (low dose, high dose). Let further n* < n, denote the 
PP sample size for group /. Finally, let T,- denote the test statistic for i — 1, . . ., 4, and p(r,-,r,) the 
correlation between T, and 7). With this notation, 

P(T.,T2) = p(T3,r.).f^)"Y^y'' 

which reduces to 0.5 if «o = «i = «2- Further, 

V «o«i J Wo+nV V «o«2 J V«o+«y 
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which both reduce to (wq/hq)''^^ for no = rii and nl — n*, i=l,2. Finally, p(ri,r4) = 
p(Tu Ti)p(T3, 74) and p(r2, ^3) = piTi, T4)p(T3, T4), which both reduce to l/2(«*/no)'/^ for no = 
n\ —ni and «q — n* — n\. In this simplest case of equal group sample sizes within PP and ITT we 
thus have, assuming Ko/«o — 0-9 as an example 

( 1 for /= {1,3}, /= {2,4} and/= {/}, /= 1,...,4 

cj = <^ 1.0783 for / = {3,4} and for all / ^ {1,2} 
y 1.0706 otherwise 

As a consequence, the resulting family of tests is no longer consonant, although the 
differences in the resulting local significance levels are small. For example, 
qi,2,3,4)Wi({l,2,3,4})a = 0.0135>0.0134 = C(i,4)H'i({l,4})a, violating condition (4). Similar to 
Example 2, we can enforce consonance by applying the graphical test procedure from Figure 4 with 
5 = 0.0071. 

Finally, we note that this multiple test procedure is immediately applicable to testing for a 
treatment effect at two different dose levels in an overall population and, if at least one dose is 
significant, continue testing in a pre-specified subpopulation. This could apply to testing, for ex- 
ample, in the global study population and a regional subpopulation or in the enrolled full popu- 
lation and a targeted genetic subpopulation. 



3.3 Weighted Simes tests 

Generalization of the original Bonferroni-based graphs from Section 3.1 also apply when the 
correlations between the test statistics are not exactly known, but certain restriction on them are 
assumed. A typical case in practice is to assume (or show) that the test statistics have a joint 
multivariate normal distribution with non-negative correlations. In this case, the Simes test is a 
popular test. Here, we discuss the use of a weighted version of the Simes test for the intersection 
hypotheses Hj,J c /. 

The unweighted Simes test, as originally proposed by Simes (1986), rejects Hj if there exists a /e/ 
such that p{j) <j/ma., where p(\) < ■ ■■ < P(m) denote the ordered /^-values for the hypotheses 
Hi, i e I. The Type 1 error rate is exactly a if the test statistics are independent and it is bounded by 
ot if positive regression dependence holds. This follows from Benjamini and Yekutieli (2001), who 
showed false discovery rate control for a related step-up procedure under positive regression 
dependence on the test statistics. Note that this condition is not always easy to verify or even justify 
in practice. 

The weighted Simes test introduced by Benjamini and Hochberg (1997) rejects if/ if for some jel 
P(j) — S=i ^(i)' where = vv^a and denotes the weight associated with /?(,). An equivalent 
condition is to reject Hj if for some jel 

Pj <'^<Xi = <x'^Wi (5) 

i€lj ielj 

where Ij — {k e /; pk < Pj}- This weighted Simes test reduces to the original (unweighted) Simes test 
if Wi = l/m,i € I. Kling (2005) showed that the weighted test is conservative if the univariate test 
statistics are positive regression dependent for any number of hypotheses. This, for example, is the 
case if the test statistics follow a multivariate normal distribution with non-negative correlations 
and the tests are one-sided (Benjamini and Heller, 2007). 

For given weights Wf(J), J^I, and assuming positive regression dependence among the univariate 
test statistics for all m hypotheses //,-, iel, the weighted Simes test can be applied to all intersection 
hypotheses Hj, J^I. By means of the closure principle the resulting multiple test procedure rejects 
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Hi, iel, at level a if for each with ieJ, there exists an index jeJ such that 

Pj<aJ2wkiJ) (6) 

keJj 

where Jj — {k e J;pk < pj}. This follows from the application of condition (5) to all subsets /c/, 
and the fact that any subset of m positive regression dependent test statistics is also positive 
regression dependent. Related gatekeeping procedures based on the Simes tests were described, 
among others, by Dmitrienko et al. (2003) and Chen et al. (2005). 

If all weights are equal, the above procedure reduces to the procedure by Hommel (1988), which 
is known not to be consonant. In case of unequal weights, a corresponding sequentially rejective test 
procedure is not available and one may have to go through the entire closed test procedure using 
weighted Simes tests for each intersection hypotheses. Nevertheless, for a given weighting strategy, 
this procedure is uniformly more powerful than an associated Bonferroni-based procedure from 
Section 3.1. This follows from the fact that any hypothesis rejected by the closed weighted Bon- 
ferroni test procedure can also be rejected by the corresponding closed weighted Simes test pro- 
cedure; see, for example, the Appendix in Maurer et al. (2011). 

Although full consonance is generally not available for Simes-based closed test procedures, we 
can still derive a partially sequentially rejective test procedure which leads to the same test decision 
as the closed test procedure defined in (6). In the following, we assume that the weights are ex- 
haustive, i.e. J2k^jWkiJ) = 1 for all subsets J el. 

Algorithm 4 (Weighted Simes Test) 

(i) If Pi>(x for all iel, stop and retain all m hypotheses. 

(ii) If pi<ci for all iel, stop and reject all hypotheses. 

(iii) Perform the Bonferroni-based graphical test procedure from Section 3.1. Let denote the 
index set of rejected hypotheses and its complement in /. If < 3, stop and retain the 

remaining hypotheses. 

(iv) If 1 7,5 1 > 3 consider the weights Wi(/^'), / e I^, and the transition matrix G defined on as 
the new initial graph for the remaining hypotheses. Compute the weights Wk(J) for all 
J c with Algorithm 1 . 

(v) Reject Hj, i e /,f, if for each / c with ieJ, there exists an index jeJ such that 

Pj<aJ2MJ)- (7) 

k€jj 



With step (ii), all hypotheses Hj, iel can be rejected if pj<a- for all /e/. This follows from 
the fact that for each / there is always a largest p,, jeJ, such that ]j = J and therefore 
^^kei ^^'k(J) — '^Y^ke/ ^^'k(J) — ^- Hence condition (6) holds for all J^I and therefore for all H,, 
iel. Note that if the weights are not exhaustive, step (ii) may no longer be valid and should be 
skipped. 

The stopping condition in step (iii), |/^| <3, is explained as follows. Assume first that | — 1, i.e. 
one hypothesis is left, say i/,. If /?, < a, one would have rejected already all hypotheses in step (ii) and 
stopped the procedure because for all other hypotheses than H, necessarily pj<i. Therefore, />,>a 
and one cannot reject H,. Similarly, if |//| = 2, the respective /^-values cannot be both smaller than a. 
Also if only one of them, say pi, is smaller and the other is larger than a, then pi > Wi(I^)a, since 
otherwise the Bonferroni test in step (iii) would have rejected //,. In that case the Simes test cannot 
reject //,• either and hence both remaining hypotheses must be retained. 

Algorithm 4 is essentially looking first for outcomes that are easy to verify (steps (i) and (ii)) or 
where sequential rejection of the hypotheses is possible (step (iii)). Only then one needs to compute 
for all remaining hypotheses and their subsets the weights and apply the closed weighted Simes 
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procedure as given in (6). It can happen though that no hypotheses can be rejected in the first three 
steps and that one has to perform step (iv) with the full set of all m hypotheses. Note that one could, 
of course, start immediately with step (iv) on the full hypotheses set. The resulting decisions 
are identical to those obtained with Algorithm 4, because for any given weighting strategy, any 
hypothesis rejected by the closed weighted Bonferroni test procedure is also rejected by the asso- 
ciated closed weighted Simes test procedure. 

Similar to the case that knowledge about the joint distribution of the p-values is partially missing 
(as discussed in Section 3.2), we consider now the case that positive regression dependence cannot 
be assumed between all m test statistics. Let 4, A = 1, ...,/< m, be a a partition (i.e., / = uj^^j//, and 

n /, = 0 for h ^ 0 such that for each family of hypotheses i/,, / e h, positive regression depen- 
dence between the respective test statistics holds. Then we can reject Hj, / c /, if for some j and h 
with j e //, — 



where Jhj — {k e Jh',Pk < Pj}- This procedure controls the Type I error rate at level a for any in- 
tersection hypothesis Hj. This is seen as follows. The weighted Simes test is applied separately to each 
of the partition sets Jh of /. With the definitions for //, and Hhj above, for a fixed e {1, ...,/), the 
probability of the event that there exists a jeJ/, such that pj < 7. J^ke/, "VfC-^). is less than or equal to 
^ YlikeJH ^k{J) by the weighted Simes test. Hence the probability that this happens in any of the 
partitions is less than Yl\i=\ '^Ylke.i,, '^'^'kiJ) = '^J2k€j^k(J) = a by means of the Bonferroni in- 
equality. For a given partition //,, h — {,..., I, with "local" regression dependence within the disjunct 
subsets of associated test statistics, condition (7) in the algorithm hence can be replaced by (8). 

We conclude this section with an example. For the weighting strategy from Example 1, the 
resulting closed weighted Simes test will reject more hypotheses than the related closed weighted 
Bonferroni test only if all four ;7-values are less than or equal to a (Maurer et al., 201 1). The latter is 
not the case for the numerical example in Section 3.1, because, for example, 773 = 0.1 > 0.025 = a and 
hence no further hypothesis can be rejected. However, if we had instead, for example, p2, = 0.015 
and /74 = 0.022, the closed weighted Simes test would reject all four hypotheses, two more than with 
the closed weighted Bonferroni test. Generally speaking, the weighted Simes test has power ad- 
vantages over alternative weighted test procedures if the effect sizes are of similar magnitude. 

4 gMCP package in R 

The gMCP package (Rohmeyer and Klinglmueller, 2011) in R (R Development Core Team, 2011) 
currently implements the Bonferroni-based graphical approach from Section 3.1 and the closed 
weighted parametric tests from Section 3.2. R is a language and environment for statistical 
computing and graphics (Ihaka and Gentleman, 1996). It provides a wide variety of statistical 
and graphical techniques, and is highly extensible. The latest version of gMCP is available 
at the Comprehensive R Archive Network (CRAN) and can be accessed from http://cran. 
r-project . org/pack:age=gMCP/. In the following, we give only a brief illustration of the gMCP 
package. We refer to the installation instructions at http://cran.r-project.org/web/ 
packages/gMCP/INSTALL and the accompanying vignette for a description of the full func- 
tionahty (Rohmeyer and Klinglmueller, 2011). 

4.1 Weighted Bonferroni tests witli gMCP 

We consider the cardiovascular clinical trial example from Dmitrienko and Tamhane (2009) to 
illustrate the implementation of the Bonferroni-based graphical approach from Section 3.1 in 
the gMCP package. The trial compared a new compound with placebo for two primary and two 




(8) 
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secondary endpoints. Consequently, we have two families of hypotheses J^\={H\,H2} and 
J^2 = {Hs,H4}. 

Dmitrienko and Tamhane (2009) used this example to illustrate the truncated Holm procedure 
described in Dmitrienko et al. (2008) and Strassburger and Bretz (2008). Given multiple families of 
hypotheses in a pre-specified hierarchical order, the key idea of truncated tests is to avoid propagating 
the complete significance level within a family until all its hypotheses are rejected in order to proceed 
testing the next family in the hierarchy. Instead, once at least one hypothesis is rejected in a given 
family, a fraction of the significance level is reserved to test subsequent families of hypotheses. In 
principle, truncation can be applied to any of the test procedures discussed in Section 3. 

In the cardiovascular study example, the hypotheses in are only tested, if at least one 
of the hypotheses in are rejected. We assume that J^i is tested using the truncated Holm procedure 
with truncation parameter ye [0,1]. Let P{i)<P{2) denote the ordered ;?-values with associated hypotheses 
//(I) and Hf2). Consequently, H(\) is tested at level a/2. If H(\) is rejected, H(2) is tested at level 
a/2+y(a/2). The family J^2 is then tested with the regular Holm procedure either at level (l-Y)a/2 or at 
level a, depending on whether only one or both hypotheses in are rejected, respectively. 

The gMCP package offers a GUI to conveniently create and perform Bonferroni-based graphical 
test procedures, such as the one for the test procedure above. To this end, we invoke in R the gMCP 
package and subsequently call the GUI with 

> library (gMCP) 
> graphGUI ( ) 

Different buttons are available in the icon panel of the GUI to create a new graph. The main 
functionaHty includes the possibiHty of adding new nodes as well as new edges connecting any two 
selected nodes. In many cases, the edges will have to be dragged manually in order to improve the 
readability of the graphs. The associated labels, weights, and significant levels can be edited directly in 
the graph. Alternatively, the numerical information can be entered into the transition matrix and other 
fields on the right-hand side of the GUI. Figure 5 displays the complete test procedure for the cardi- 
ovascular study example using the gMCP package: The truncated Holm procedure for T\ with trun- 
cation parameter y and the regular Holm procedure for 7^2- Note that we can immediately improve that 
test procedure by connecting the secondary hypotheses and H4 with the primary hypotheses H\ and 
H2 through the s-edges introduced in Bretz et al. (2009). We refer to the vignette of the gMCP package 
for a description of how to construct 8-edges with the GUI (Rohmeyer and Klinglmueller, 2011). 

The GUI offers the possibility to perform sequentially Bonferroni-based test procedures defined 
through a graph hke the one displayed in Figure 5 and in addition to calculate adjusted /(-values as 
well as simultaneous confidence intervals. To illustrate this functionaUty, we consider Scenario 1 
from Dmitrienko and Tamhane (2009) and assume the unadjusted /^-values px = 0.0121, = 0.0337, 
;?3 = 0.0084, and /74 = 0.0160, which are entered directly into the GUI. By clicking on the corre- 
sponding button in the icon panel and and specifying y = 0.5, one obtains in this example the 
adjusted /7-values 0.024, 0.045, 0.045, and 0.045 for the four hypotheses Hy, H2, H^, and H4, 
respectively. These adjusted /i-values are identical to those reported in Dmitrienko and Tamhane 
(2009). Accordingly, one can reject all four hypotheses at level a = 0.05. Simultaneous confidence 
intervals can be obtained as well from the GUI after entering additional information on effect 
estimates and standard errors. Finally, the user may perform the sequential test procedure by 
clicking on the green triangle in the icon bar. By doing so, the "Reject" buttons in the lower right 
become activated and one can step through the graph as long as significances occur. 



4.2 Weighted parametric tests with gMCP 

The gMCP package provides also a convenient interface to perform graphical test procedures 
without the GUI using the R command fine. We illustrate this with the closed weighted parametric 
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gMCP GUI 



File Example graphs Analysis Extras Help 



pv4 



Place new nodes and edges or start the test procedures 



Description Analysis 




Bonferroni-based test procedure for the cardiovascular trial example from 

□ mitrienko and Tamhane (2009). 



Sum of weights: 1 

Total a : 0.05 



• No Information about correlations 
Select a standard correlation 
I Select an R correlation matrix 



1 No quadratic matrices found. 



1" ^ 


^^^J j 0.0121 


1 Reject and pass ct 1 






1 I 0.0337 


[ Reject and pass a 1 




1° 


t 1 0.0084 


1 Reject and pass a 1 




1° 


I 1 0.0160 


I 1 Reject and pass a | 



Figure 5 Screensliot of the GUI from the gMCP package. Left: Display of the graphical 
Bonferroni-based test procedure for the cardiovascular trial. Right: Transition matrix, 
initial weights and unadjusted values. 



tests from Section 3.2 and revisit Example 2. We first define the related transition matrix G and the 
weights H',(/), iel, through 

> G <- matrix (0, nr=4, nc = 4 ) 
>G[1,3] <-G[2,4] <-G[3,2] <-G[4,l] <-l 
>w <- c(l/2, 1/2, 0, 0) 

The function matrix2graph then converts the matrix G and the vector w into an object of type 

graphMCP 

> graph <- matrix2graph (G, w) 

> graph 

A graphMCP graph 

Overall alpha: 1 

Hi (not re j ected, alpha=0.5) 

H2 (not re j ected, alpha=0.5) 

H3 (not re j ected, alpha=0) 

H4 (not re j ected, alpha=0) 

Edges : 

HI - (1) -> H3 

H2 - (1) -> H4 

H3 - (1) -> H2 

H4 - (1) -> HI 
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The gMCP function takes objects of the type graphMCP as its input together with a vector of 
j7-values and performs the specified multiple test procedure. In particular, one can specify a cor- 
relation matrix with the effect that a closed weighted parametric multiple test procedure is per- 
formed under the standard analysis-of-variance assumptions with known common variance. 

In Example 2 we assumed normally distributed test statistics with a block-diagonal correlation 
matrix of the form 

/ 1 0.5 NA NA\ 

0.5 1 NA NA 

NA NA 1 0.5 ' 

\NA NA 0.5 1 / 

where NA reflects the fact that the correlation between the primary and secondary endpoints is 
unknown. Accordingly, we let 

> cr <- matrix (NA, nr = 4, nc = 4) 

> diag (cr ) < - 1 

>cr[l,2] <-cr[2,l] <-cr[3,4] <-cr[4,3] <-l/2 
and define the unadjusted /j-values 

> p <- c(0.0131, 0.1, 0.012, 0.01) 

Finally, we perform the closed weighted parametric test at a specified significance level a = 0.025, 
say, by calhng 

> res <- gMCP(graph, p, corr = cr, alpha = 0.025) 

This returns an object of class gMCPResult providing information on which hypotheses are re- 
jected 

> resSrejected 

Hi H2 H3 H4 
TRUE FALSE TRUE FALSE 

We conclude from the output that both Hx and Hj, can be rejected. We come to the same con- 
clusions, if we report the adjusted ;7-values and compare them with a = 0.025 

> resSadj PValues 

HI H2 H3 H4 

0.02431856 0.10000000 0.02431856 0.10000000 

Alternatively, one can use a sequentially rejective Bonferroni-based test procedure from Section 3.2 
by omitting the corr argument 

> gMCP (graph, p, alpha = 0 . 02 5) ©rejected 

HI H2 H3 H4 
FALSE FALSE FALSE FALSE 

As seen from the output, none of the null hypotheses can be rejected, which coincides with our 
conclusions from Section 3.2. 
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5 Discussion 

This paper shows that the graphical approach introduced by Bretz et al. (2009) can be used to create 
and visualize tailored strategies for common multiple test problems. By dissociating the underlying 
weighting strategy from the employed test procedure, it is seen that the graphical approach is not 
restricted to Bonferroni-based tests. Similarly, the graphs introduced by Burman et al. (2009) define 
weights for ah intersection hypotheses and the procedures discussed in this paper can be apphed 
using these weights. Extended graphical approaches include weighted Simes tests and weighted min- 
p tests in the sense of Westfall and Young (1993). The latter take into account all or some of the 
joint multivariate distributions of /7-values. Consonance and the corresponding shortcuts may be 
lost, but for any concrete multiple test strategy, consonance can be checked prior to a clinical study. 
As shown in this paper, consonance can be enforced and related sequentially rejective graphs 
established at least in some simple situations. 

Many proposed multiple test procedures in the literature can be expressed with the methods 
described in this paper. On the other hand, the methods in this paper also allow one to investigate 
alternative procedures that go beyond the published results. But even if the closure principle is very 
common in practice, it does not necessarily lead to consonant multiple test procedures. We gave 
monotonicity conditions for ensuring consonant graphical weighting strategies, but it is not always 
clear when these conditions are met if weighted parametric or Simes tests are used. In principle, one 
could enforce consonance following, for example, the approach of Romano et al. (2011), although 
the computation of the rejection regions could become tedious. We leave this topic for further 
research. 
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